The 1.4 GHz light curve of GRB 970508 
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ABSTRACT 

We report on Westerbork 1.4 GHz radio observations of the radio counterpart 
to 7-ray burst GRB 970508, between 0.80 and 138 days after this event. The 1.4 
GHz light curve shows a transition from optically thick to thin emission between 
39 and 54 days after the event. We derive the slope p of the spectrum of injected 
electrons (dA/d7 e oc % p ) in two independent ways which yield values very close 
to p = 2.2. This is in agreement with a relativistic dynamically near-adiabatic 
blast wave model whose emission is dominated by synchrotron radiation and in 
which a significant fraction of the electrons cool fast. 

Subject headings: gamma rays: bursts — gamma rays: individual (GRB 970508) 
— radio continuum: general 

1. Introduction 

The peak luminosities of 7-ray bursts (GRBs) are highly super-Eddington and require 
relativistic outflows (Paczyhski 1986; Goodman 1986). Paczyhski and Rhoads (1993) 
pointed out that radio emission is expected as a result of the interaction between such 
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a relativistic outflow and an external medium, as is, e.g., observed in extragalactic jet 
sources (see also Katz 1994; Meszaros and Rees 1997). They estimated that the strongest 
GRBs may be followed by transient (~ 10 mJy) radio emission at intervals ranging from 
minutes (for a distance d ~ 10 5 pc) to several weeks (d ~ 10 9 pc). However, searches 
for radio counterparts through follow-up observations (Frail et al. 1994,1997a; Koranyi 
et al. 1995; Galama et al. 1997a,b) were without success until recently. With the rapid 
accurate location capability of the Wide Field Cameras (WFCs; Jager et al. 1995) onboard 
the Italian-Dutch X-ray observatory BeppoSAX (Piro et al. 1995) it has recently become 
possible to detect fading X-ray, optical and radio counterparts to GRBs (Costa et al. 1997a; 
Piro et al. 1997a; Groot et al. 1997a,b; Van Paradijs et al. 1997; Galama et al. 1997c,1998a; 
Sahu et al. 1997; Bond 1997; Metzger et al. 1997; Frail et al. 1997b,c; Bremer et al. 1998; 
Halpern et al. 1997). These observations have settled the discussion on the GRB distance 
scale ('galactic halo' versus 'cosmological', see e.g. Fishman and Meegan 1995, Lamb 1995, 
Paczyhski 1995): GRBs occur at Gpc distances. 

GRB 970508 is the first GRB to be detected in the radio (Frail et al. 1997b,c); the 
radio source position coincides with that of the optical (Bond 1997) and X-ray (Piro et 
al. 1997b) afterglow sources. Assuming that the variations of the source at 4.86 and 8.46 
GHz are due to interstellar scintillation (ISS), their damping with time is consistent with 
a highly relativistically expanding shell passing a diameter of ~ 3/ias (Frail et al. 1997b). 
VLBI observations show that the source is unresolved (< 0.3 mas, Taylor et al. 1997). 

We here report on the results of 1.4 GHz radio observations of GRB 970508, made 
with the Westerbork Synthesis Radio Telescope (WSRT) between 0.80 and 138 days after 
the burst occurred. 



2. Radio Observations 

On May 8.904 UT BeppoSAX recorded a moderately bright GRB (Costa et al. 
1997b) with the Gamma-Ray Burst Monitor (GRBM; Frontera et al. 1991), which was 
also recorded with the Wide Field Cameras on board BeppoSAX. Analysis of the WFC 
observations gave a 3' (3a radius) error box centered on RA= 06 h 53 m 28 s , Dec = +79°17f4 
(J2000; Heise et al. 1997b). The burst was also recorded (Kouveliotou et al. 1997) with 
the Burst and Transient Source Experiment (BATSE; Fishman et al. 1989) on board the 
Compton Gamma-Ray Observatory. 

The error box of GRB 970508 was first observed at 1.4 GHz with the WSRT at the 
preliminary position given by Heise et al. (1997) on May 9.70 UT, starting 0.80 days 
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after the event, for 9.0 hours. We used the standard WSRT 1.4 GHz receiving system and 
continuum correlator, providing us with 5 bands of width 10 MHz and 3 bands of width 5 
MHz. The noise level in a continuum map after 12 hours of integration is typically 0.05 
mJy beam -1 for 14 telescopes and the full 65 MHz bandwidth. For the declination of GRB 
970508, at 1.4 GHz, the synthesized beamwidth is 13" x 13" (full width at half maximum; 
FWHM), and the field of view is about 0?6. Table [l] provides a log of the observations. 

The data were analyzed using the NEWSTAR software package^. The interferometer 
complex visibilities were first examined for possible electromagnetic interference and other 
obvious defects. Bands with strong interference were either deleted or carefully edited. 
Interference was usually limited to about 5 % of the data. The interferometer gain and 
phase, for each observation and each band, were calibrated using the standard WSRT 
calibrators 3C 48, 3C 147 and 3C 286 (15.96, 9.50 and 14.77 Jy, respectively at 1.4 GHz; 
Baars et al. 1977). We constructed a model of the GRB 970508 field (used for deconvolution) 
that contains 95 point sources and 60 clean components, (together about 100 sources) and 
is complete down to a level of ~ 0.45 mJy. No self calibration was performed. 

2.1. The light curve 

Due to noise the location of a source in the radio map can shift by roughly 0.5 
FWHM/SNR, where SNR is the signal-to- noise ratio of the detection. For very low SNR 
detections the shift can be fairly large. If we would take a peak in the map near the location 
of GRB 970508 we would bias ourselves systematically to higher flux densities (as we cannot 
discern whether the peak is due to noise on top of the detection or the detection itself). 
Therefore, we constructed the light curve using the flux density at the pixel on the exact 
position of the radio counterpart (Frail 1997c). The errors (la) are the r.m.s. map noises 
(determined by a quadratic fit to the deconvolved image). The maps are super sampled 
with ~ 6 pixels beam -1 . Negative flux density values in the 1.4 GHz light curve are due to 
r.m.s. noise fluctuations on a non- or weakly-detectable source. The 1.4 GHz light curve is 
shown in Fig. [l|. We also divided the data from May 9 until September 23 into 8 parts and 
added the observations in each part to determine average flux densities. In Table ^| and in 
Fig. [I] we present these average flux densities. 

During the first 40 days after the onset of GRB 970508 the 1.4 GHz radio counterpart 
is not detected. A combined map of all observations (May 9. 7- June 16.44 UT; excluding 
June 8.43 UT for reasons of data quality) yields 33 ± 40 yuJy. The average 4.86 GHz 
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flux density during the first 50 days is 560 ± 10 /xJy (Frail et al. 1997b), implying that, 
on average, the spectral index ai.4-4.86GHz > 1-5 (where we used a 2a upper limit of 80 
/ijy at 1.4 GHz), i.e. the radiation is self-absorbed (F v cx z/ 2 ; Katz and Piran 1997). At 
later times (t > 50 days) 1.4 GHz emission is detected at, on average, 228 ± 30 /iJy (Jul 
1.49-Sep 23.22; 53.6-138.3 days after the event). From Tab. § we see that then the spectral 
indices ai.4_4.86GHz and cm.4-8.46GHz are consistent with the expected low frequency tail of 
synchroton radiation (F u oc z/ 1 / 3 ; Rybicky and Lightman 1979). Hence a transition from 
optically thick to thin emission occurred between 39 and 54 days after the event (see also 
Frail et al. 1997b). 

During the first month the 8.46 and 4.86 GHz flux densities show rapid fluctuations, 
attributed to ISS, with an ISS diffractive time scale t<jif of less than a day, and decorrelation 
bandwidth z/d ec of less than 2 GHz (Frail et al. 1997b). For 1.4 GHz this implies tdif ~ 
3 hours and z/d ec ~ 0.8 MHz (Goodman 1997). Hence ISS is not likely to affect the 1.4 
GHz light curve: it was constructed from observations with long integration times (> 3 
hrs) and the full bandwidth (65 MHz). At 1.4 GHz the critical angular size below which 
diffractive scintillation can be observed is a factor 4.5 smaller than at 4.86 GHz (Goodman 
1997). Assuming that the apparent expansion velocity of the blast wave is constant and 
using the fact that Frail et al. (1997b) observed scintillation at 4.86 GHz until day ~ 30 
we would expect diffractive ISS to modulate the 1.4 GHz flux density for t < 7 days. We 
have searched for diffractive ISS in the individual 5 and 10 MHz bands during the first two 
weeks, dividing the data into segments of 6 hours (a trade off between the ISS decorrelation 
time and a sufficient amount of data). We find only one > 3a detection (550 ± 160 /iJy; 3.4 
a) on May 11.51 at 1.415 MHz (10 MHz bandwidth). The probability of a 3.4 a detection 
in ~ 100 independent measurements (with an asumed Gaussian distribution) is ~ 3 %. We 
conclude that ISS was not observed by us at 1.4 GHz. 

3. Discussion 

The observed optical spectral slope a and the optical power law decay of the light 
curve F v oc t s is not consistent with the expected relation for the simplest blast wave model 
(5 = 3a/2; e.g. Wijers, Rees and Meszaros 1997). The observed power law decay value, 
5 = -1.141 ± 0.014 (t > 2 days, Galama et al. 1998b; see also Pedersen et al. 1998, 
Castro-Tirado et al. 1998, Sokolov et al. 1998) would imply a = —0.761 ± 0.009, while 
in the optical passband a = —1.12 ± 0.04 is observed (Galama et al. 1998c, from here on 
Paper II; see also Sokolov et al. 1998). In the following we show that this may be explained 
by rapid cooling of a significant fraction of the electrons. 



- 5 - 



A population of electrons with a power-law distribution of Lorentz factors 7 e 
(cL/V/oV/c oc 7~ p ) above some minimum value 7 m emits a power law synchrotron spectrum 
above the frequency u m (corresponding to radiating electrons with 7 m ; e.g. Rybicki & 
Lightman 1979). Independently, above some Lorentz factor 7 C the electrons may cool 
rapidly, and an extra break in the spectrum is expected at the corresponding frequency v c 
(Sari, Piran, & Narayan 1998). Beyond a certain time t (t is small ~ 500 sec; see Paper 
II) the evolution of the blast wave is adiabatic (Sari et al. 1998 and see Paper II); then 
z/ m < v c , and the spectrum varies as F v oc u^^ 1 ^ 2 from u m up to z/ c ; above u c it follows 
F u oc v~ v l 2 and below u m it follows the low frequency tail, F v oc z/ 1 / 3 (Sari et al. 1998). 
The evolution in time of the GRB afterglow is determined by the evolution of these break 
frequencies: u c oc t~ 1 / 2 and u m oc t~ 3 / 2 (both decrease with time). 

The decay part of the optical R c (Coussins R) band light curve (in the optical passband 
v > v c for t ^ 1.2 days; Paper II) goes as F v oc t( 2_3 P)/ 4 , while the spectrum is then 
F v oc v~ v l 2 (Sari et al. 1998). This allows us to make two independent measurements 
of p: using F Rc oc r 1 - 141 * - 014 we find p = 2.188 ± 0.019 and using a opt = -1.12 ± 0.04 
gives p = 2.24 ± 0.08. The excellent agreement between the values of p supports that a 
significant fraction of the electrons cool rapidly and that the evolution of the GRB remnant 
is adiabatic. Additional evidence for rapid cooling of a significant fraction of the electrons 
is given in Paper II. 

Observations by Bremer et al. (1998) with the IRAM Plateau de Bure Interferometer 
(PdBI) at 86 GHz show a maximum around ~ 12 days. We identify this maximum with the 
break frequency v m passing 86 GHz at t m ,86GHz ~ 12 days (Paper II). We expect, the 8.46 
and 1.4 GHz emission to peak at t mj 8.46GHz ~ 55 days and t mj i.4GHz ~ 180 days, respectively 
(z/ m oc t~ 3 / 2 ). Near day 55, a shallow maximum can be seen in the 8.46 GHz light curve 
(Frail et al. 1997b). Unfortunately our 1.4 GHz light curve cannot be used to test the 
presence of the maximum at that frequency, both due to low signal to noise and because it 
ends 150 days after the burst, i.e. before the predicted maximum. 

Before z/ m passes 8.46 GHz at t mj 8.46GHz we expect the 8.46 GHz spectrum to follow the 
low frequency tail F v oc z/ 1 / 3 , while after t m) 8.46GHz it is expected to be F v oc u^^ 1 ^ 2 = z/~ - 6 
(Sari et al. 1998 and we have used p = 2.2). Thus, we predict a gradual transition between 
£ m ,8.46GHz ~ 55 days and t m ,4.86GHz ~ 80 days (when also at 4.86 GHz z/ m has passed) from a 
= 1/3 to a = —0.6. We note that this expectation is different from blast wave models that 
do not include the effect of rapid cooling of a significant fraction of the electrons (F u oc v~ 1A 
similar to the optical slope; see e.g. Wijers et al. 1997). Also the decays at 8.46 GHz (after 
i m ,8.46GHz ~ 55 days) and 4.86 GHz (after t m ,4.86GHz ~ 80 days) are expected to be different 
from that in the optical and X-ray passbands, F v oc t 3(1 ~p)/ 4 = t~ - 9 ; where we have used p 
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= 2.2). These predictions can be tested with the continued monitoring at the VLA at 4.86 
and 8.46 GHz by Frail et al. (1998). 

The radio afterglow light curves of GRB 970508 (Frail et al. 1997b and this Letter) 
show a much more gradual evolution than expected (see e.g. the fit to the 8.46 and 
4.86 GHz data by Waxman, Kulkarni and Frail 1998). Also a constant self-absorption 
frequency was expected (e.g. Waxman et al. 1997) while we here show that a transition 
from optically thick to thin emission occurred around ~ 45 days. For t < t Sari et al. 
(1998) predict a decrease with time of the self-absorption frequency, u a , while for t > to the 
self-aborption frequency u a remains constant. The transition from optically thick to thin 
1.4 GHz radiation then suggests that t ~ 45 days. Also the 8.46 GHz light curve (Frail et 
al. 1997b) suggests that t cannot be much smaller than 10 days, i.e. 10 days < t ^ 55 
days (we have extrapolated backwards in time from the 8.46 GHz peak at t m ~ 55 days 
with the expected dependence F u oc t 1 / 2 for times t < t m ). This is not in agreement with 
the finding that to ~ 500 sec (Paper II). However, the absence of a break in the smooth 
power law decay of the optical light curve from 2 to 60 days after the burst (Pedersen 
et al. 1998; Castro-Tirado et al. 1998; Sokolov et al. 1998; Galama et al. 1998c) shows 
that there is no important transition in that period. This does imply that some additional 
ingredient is needed; for example, Waxman et al. (1998) argue that the transition from 
ultrarelativistic to mildly relativistic expansion of the blast wave may explain the decrease 
in the self-absorption frequency z/ a with time and the slow time dependence of the early 
radio light curves. 

The excellent agreement in the derived value for p (p = 2.2) from the decay of the 
optical light curve and the optical spectral slope support an adiabatic dynamical evolution 
of the GRB remnant and an extra break in the synchrotron spectrum at the frequency u c 
above which the radiation is from electrons which cool rapidly compared to the remnant's 
expansion time. We predict a transition in the radio spectral index 04.86-8. 46GHz from 1/3 
to -0.6, between 55 and 80 days; the light curves are predicted to decay as F u oc t~°' 9 after 
55 days at 8.46 GHz and 80 days at 4.86 GHz. 
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Table 1: 1.4 GHz WSRT observations. 
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Table 2: 1.4, 4.86 and 8.46 GHz long-term averages and 1.4 GHz spectral indices (with 
respect to 8.46 and 4.86 GHz). The 8.46 and 4.86 GHz averages have been determined from 



observations by Frail et al. (1997b). 
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1,4 GHz lightcurve of GRB970508 








t (days after burst) 



Fig. 1.— The 1.4 GHz light curve of GRB 970508. The data are from Tab. | (•) and Frail et 
al. (1997b; * and the 2a upper limits). The histogram represents long-term average WSRT 
1.4 GHz flux densities from Tab. §. 



